Probability distribution of the vacuum energy density 
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As the vacuum state of a quantum field is not an eigenstate of the Hamiltonian density, the vacuum 
energy density can be represented as a random variable. We present an analytical calculation of 
the probability distribution of the vacuum energy density for real and complex massless scalar fields 
in Minkowski space. The obtained probability distributions are broad and the vacuum expectation 
value of the Hamiltonian density is not fully representative of the vacuum energy density. 



PACS numbers: 11.10.-z; ll.10.Ef; 03.70.+k; 98.80.Cq 



O 

(N 

> 
O 

(N 



(N 
> 

00 

c4 
o 
o 



X 



Quantum field theory (QFT) is a theoretical frame- 
work that has enabled, via its phenomenologically suc- 
cessful models such as the Standard Model, a detailed 
quantitative understanding of a vast number of physical 
phenomena. In numerous areas of its application QFT 
predictions are confirmed at ever greater accuracy with 
every new experiment. Yet, these QFT models with such 
a remarkable performance in precise determination of 
scattering cross-sections and decay rates also give a defi- 
nite prediction for the vacuum expectation value (VEV) 
of the Hamiltonian density. As long as we disregard grav- 
itational interaction, the size of this VEV is irrelevant [l[ . 
When the gravitational interaction, described by general 
relativity, is included the VEV of the Hamiltonian den- 
sity can no longer be ignored. Moreover, any estimate of 
this quantity within QFT, also referred to as zero point 
energy density, shows that this contribution to the source 
term of the Einstein equation is crucial. The zero point 
energy density represents a major contribution to the to- 
tal value of the cosmological constant (CC) @. The fact 
that the size of this contribution, as well as the size of 
other identified contributions to the cosmological con- 
stant, are many orders of magnitude larger than the ob- 
served value of the CC is at the core of the notorious old 
cosmological constant problem. A deeper understanding 
of all contributions to the cosmological constant is clearly 
needed to resolve the old CC problem, perceived as one of 
the most serious problems in theoretical physics. Another 
important physical phenomenon where the understand- 
ing of the energy density in the vacuum is paramount is 
inflation 0]. 

In this paper our primary goal is understanding of the 
energy density that quantum fields develop in the vacuum 
state. Hereafter we call this energy density the vacuum 
energy density and investigate how well it is represented 
by a VEV of the Hamiltonian density. The central result 
of this paper is the calculation of the probability density 
function (p.d.f.) of the vacuum energy density for both 



'Electronic address: gorand@thphys.irb.hr 
^Electronic address: dglavan@dominis.phy.hr 
t Electronic address: shrvoje@thphys.irb.hr 



real and complex free massless scalar fields in Minkowski 
space. We present a fully analytically tractable treat- 
ment of the problem and finally provide the analytical 
expressions for the probability density functions of the 
vacuum energy density. As our motivation comes from 
cosmology we do not assume normal ordering during cal- 
culation. However, since normal ordering is a widely ac- 
cepted regularization procedure in QFT we present "nor- 
mally ordered" results in the discussion. Throughout the 
paper we adopt the notation and conventions of [l|. 

We start with a direct observation that, for both real 
and complex scalar fields, the Hamiltonian H and the 
Hamiltonian density % operators do not commute. The 
vacuum state |0) is not an eigenstate of T-L. The mea- 
surement of the Hamiltonian density in vacuum can re- 
sult in different values of vacuum energy density with 
different probabilities. A quick calculation of the vari- 
ance o\ = (H 2 ) - (H) 2 shows that = 2/3(H} 2 for 
real and er^ = 1/3 (T-L) 2 for complex massless scalar field. 
The result that an is of the same order as (T-L) reveals 
the existence of a broad p.d.f. for the vacuum energy 
density. At this point, one can question the relevance of 
H as an observable. Namely, from the practical point 
of view, it is more natural (i.e. physical) to observe a 
smoothed energy density in a finite space-time region. In 
that case, the relevant region depends on the size of the 
system under investigation. It can be as big as the hori- 
zon in the case of cosmology or as small as experimental 
setup in the case of Casimir effect. The upper limit cor- 
responding to the infinite region gives an uninteresting 
result (p.d.f. is a delta function because % cx H) and 
the lower limit corresponds to taking % as an observ- 
able. Smoothing of energy density in a finite space-time 
region considerably complicates the calculation. Here we 
assume that % is a legitimate observable and our result 
is directly useful as a limiting value for making measure- 
ments in smaller and smaller volumes. It is important to 
note that treating % as an observable can be crucial in 
understanding cosmological inflation where our patch of 
the Universe originates from a virtually pointlike element 
of space. 

Since the spectrum of the Hamiltonian density opera- 
tor is positively defined, the p.d.f. of the vacuum energy 
density is to be determined by using a Laplace transform. 
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If p(e) is the p.d.f. of interest, its Laplace transform is 
given by 



L{p(e)} = 



e p(e)de = 
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where (H n ) is the VEV of n-th power of the Hamilto- 
nian density operator. Dimensions of s and e are [J?] -4 
and [E] , respectively. Since the original function is 
uniquely defined by its image, once the transform L{p(e)} 
is known, the p.d.f. p(e) can be determined by using var- 
ious properties of the Laplace transform 

To obtain the transform L{p(e)} one has to calculate 
(H n ), n = 1, . . . , oo. Hamiltonian densities of real (Hn) 
and complex (He) massless scalar fields are given by 
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and aL & P (a p , 6 p ) are the usual creation (annihilation) 
operators of scalar fields. For further considerations it is 
convenient to know contractions and commutators of the 
above mentioned operators. The relevant contractions 
and commutators are 
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where, from now on, arguments (x, t) of the operators r 
and c are omitted for the sake of notation simplicity. 

From Eqs. and ©, it follows that vacuum expec- 
tation value of H n can be written as 
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where the notation for momenta is renamed for later con- 
venience. 

The expectation values in integrands of Eqs. (flUf and 
(fTTj) can be calculated by using Wick's theorem [5j and 
contractions from Eqs. © - ®. It follows that 
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where P are permutations. The total numbers of terms in 
Eqs. ([11} and (HSJ) are (2n)!/(2"n!) and n\, respectively. 

Because of the symmetry on the exchange of dummy 
variables k^, it is possible to further simplify Eqs. (fTUj) 
and (jlip using binomial coefficients 
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where we assume that in the case of an ill-defined prod- 
uct dliLi an( l EliLn+l) ^ ne wn °l e contribution in corre- 
sponding brackets is equal to 1. 

Since the factors k^ • kj are odd functions of momenta 
and taking into account Eqs. (Ti"2"j) and (fT3"|) . the nonva- 
nishing contributions in Eqs. (|T4")) and (IT51) can be fac- 
torized. In factorized form Eqs. (TT4")) and (IT51 are given 
by 
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where, as before, in the case of ill-defined integrals (cases 
I = and I = n) we assume that the contribution of the 
integral is equal to 1. 

Because of the delta functions emerging from matrix 
elements, integrals in the square brackets of Eqs. (|16[) 
and (|17p reduce to the combinations of the integrals 
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where, for further convenience, a short notation (k) is 
introduced for the first integral and the second integral 
is represented in terms of the first integral. It is im- 
portant to note that the above integrals are divergent, 
so is parameter k, and strictly speaking one has to reg- 
ularize integrals flT8"|) and (fT9|) . Clearly the physically 
relevant conclusions should not depend on the choice of 
regularization. However, it is important to discuss var- 
ious regularizations since, if blindly used, they can give 
misleading results. That is especially true for the case of 
dimensional regularization. Namely, if dimensional reg- 
ularization is applied in a usual way the parameter k 
is put to zero by convention using the argument that 
it has no scale dependence. Consequently, that puts all 
the moments (except the zeroth one) of the distribution 
to zero and the p.d.f. would equal the delta functional. 
However that would imply that the Hamiltonian den- 
sity and the Hamiltonian commute (at least when their 
commutator acts on the vacuum state) which is not the 
case, as one can check by explicit calculation. To re- 
solve this paradox, one has to remember that rigorously 
speaking the integrals (ITS)) and (IT91 are not zero in di- 
mensional regularization. They are either UV divergent 
or IR divergent depending on the concrete dimension 
used for regularization. That is because relevant inte- 
gral in dimensional regularization can be written in the 
form Jq 00 dk k a (where k is magnitude of momenta in D- 
dimensional space) which is obviously divergent for any 
a. However, for practical purposes (calculation of cross- 
sections and decay rates) these integrals are convention- 
ally put to zero because their contributions are always 
fully absorbed into renormalized quantities and we have 



complete control over them. There is not any justifica- 
tion to use such an approach to calculate the p.d.f. of 
the vacuum energy density. In a sense, the standard di- 
mensional regularization procedure is not an adequate 
regularization procedure here. However, since relevant 
conclusions should not depend on the choice of the reg- 
ularization, for our purposes, we can adopt more conve- 
nient regularization. For example, in three-dimensional 
cutoff regularization where the momentum cut-off is A, 
the parameter n is equal to A 4 /(16 7r 2 ). Anyway, let us 
note that, fortunately, to obtain the functional forms of 
the probability density functions it is not necessary to 
exactly calculate the above integrals i.e. to know the 
value of «. Namely, when normalized and expressed in 
terms of expectation value the p.d.f.s do not depend on 
k explicitly (see Eqs. (|3"rJ)l and (|3"7l) later in the text). 

After integration and counting of identical terms, it 
can be shown that VEVs {W^) and (Hq) are equal to 
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where are unsigned Stirling numbers of the first 

kind. The number \S^\ gives a number of permutations 
of i elements which have exactly j cycles Q • In the above 
equations, terms are still grouped in square brackets to 
indicate the origin of each contribution if compared to 
Eqs. (fTB ]) and (fT7f . 

Now it is possible to calculate Laplace transforms of 
the probability density functions of the vacuum energy 
for real and complex massless scalar fields. After putting 
Eqs. (|20|) and (|21|) in Eq. ([T]) and rearrangement of 
terms, it follows 
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it can be seen that functions c equal 
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Finally, by taking into account identities (|25|) . (|26j) . and 
(|22|) . the relevant Laplace transforms are 

L{p n (e)} = (l + KS /3)^ 3/2 (l + ^)- 1/2 , (27) 
L{Pc(e)} = (1 + K8/S)- 3 (1 + K8)- 1 . (28) 

By using the differentiation theorem of Laplace trans- 
form, it further follows 

L{ P R(e)} = ^(-^(i + «/3r 1/2 )(i + «r 1/2 
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where L{f(c, a)} = (1 + as) x / 2 . Since 
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where theorems of linearity and convolution are used. 
After calculating the remaining integral we obtain 
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where I n (x) is a modified Bessel function of the first kind. 
Similarly, for the complex scalar field 

97 e 1 1 
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FIG. 1: The probability density functions of the vacuum en- 
ergy density for real (solid) and complex (dashed) massless 
scalar fields. Variable x is equal to the ratio e/e where e is 
the energy density and e is the expectation value which is 
different for the real and the complex cases. Mean values for 
both curves are x = 1 and the standard deviation for the 
real case is a = ^/2/3 = 0.816 and for the complex case is 

a — v^l/3 = 0.577. For the real field maximal probability is 
at x = 0.419 and for the complex field at x = 0.664. 



and finally after integration we get 
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It is convenient to express the probabilities pr,c (e) de in 
terms of vacuum expectation values 6r ; c = (Hr,c) a s 
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As can be seen from Fig. [TJ probability density func- 
tions Pr,c(x) are broad. 

The distributions given in (|3T1) and (|3~4"f depend on 
formally divergent quantity k and a regularization fol- 
lowed by renormalization must be employed. The only 
required renormalization on free field theories within the 
standard approach in QFT is given by normal ordering. 
As normal ordering corresponds to subtraction of (H) 
from "H, the effect of this renormalization on probability 
distributions is their shift so that their expectation val- 
ues vanish as shown in Fig. [5] This can be checked by 
direct calculation if : T~L := H — (Ti) is used instead of 
H, which requires calculation of (: H : ra ) in intermediate 
steps. However, the form of probability distributions for 
: % : still depends on formally divergent k. In a sense, 
this indicates that the normal ordering is insufficient to 
remove all divergences from the distributions. However, 
it is not realistic that k is really divergent. Any plausi- 
ble energy cut-off makes it big but finite. All moments 
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FIG. 2: The probability density functions of the normally 
ordered vacuum energy density for real (solid) and complex 
(dashed) massless scalar fields. Variable x is equal to the ratio 
e/n for the real and e/2n for the complex field. Parameter e is 
the energy density and k = A 4 /(16 7r 2 ) in three-dimensional 
cutoff regularization, where the momentum cutoff is A. 

can be naively renormalizcd if moments are defined as 
(: n n :) instead of (H n ) or (: U :"). Here, it is impor- 
tant to note that such procedure corresponds to infinitely 
independent renormalizations of moments since each mo- 
ment has a different divergency (beside the fact they all 
depend on the same formally divergent integral k). With 
such an approach, all the moments (except the zeroth 
one) would be equal zero and the p.d.f. would equal the 
delta functional, which leads us into contradiction with 
the fact that the Hamiltonian density and Hamiltonian 



do not commute as mentioned above. 

In conclusion, the calculations presented in this paper 
show that the vacuum energy density of a massless scalar 
field is a random variable with a broad p.d.f. Both the 
expectation of the vacuum energy density and the most 
probable value of the vacuum energy density are not rep- 
resentative for the distribution as a whole. There exists 
a nonnegligible probability for a value of vacuum energy 
that differs a lot from e. These results, obtained analyti- 
cally and in a closed form, imply that the vacuum energy 
density is a more complex object than just the VEV of 
H. 

The central results of this work have implications be- 
yond the scope of this paper. A thorough understanding 
of the vacuum energy density is essential for the expla- 
nation of the epochs of accelerated expansion such as the 
inflationary phase or the late-time accelerated expansion. 
An especially interesting prospect is that a multiverse 
scenario could be realized in the framework of QFT. To 
this end, the present analysis needs to be generalized to 
curved space-times to study the zero point energy den- 
sity in the presence of gravity. Natural extensions of the 
present analysis comprise the study of effects of mass, 
other types of quantum fields, and field interactions. The 
analytical results presented here represent a sound basis 
for the exploration of these directions of research. 
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